^ i NASA 



LOAN COPY: RETURN TO 
AFWL TECHNICAL LIBRARY 
KIRTLAND AFB, N.M. 87117 


Technical 

Report 


distributed by 

Defense Technical Information Center 

DEFENSE LOGISTICS AGENCY 

Cameron Station • Alexandria, Virginia 22304-6145 


UNCLA 



I F I E D 


TECH LIBRARY KAFB, 





TECH LIIHABY KAFB, HM 





;3beci and tawren 




fferch 1SS1 





Nasons! AefcnaijtiD; atid 
Sctace Administration 


United Gtaiffis ,A'my T 
Avjation Rgssarctt ' ’.''ij < ■■] 
Ej'iiJ Da-'c-opmant f, 
Convsrard 





fv''’5i*5A Tec^fr.;;;?■ J^gretisrafso'nrs'’ SltSS 






■ "N. 






,Sf3'/^ V t'l? i 

I k 

'»''.{»■■'■ Ai .>':?•*' a*>':' 

■^,,^ v^jli- «-,« 

V :. 


' I 


3 © B w w TO 'W' 





'i|53 


tf^.CT-’.aWWW’gOBK WW ’****'* 


Tuncsr Csbsci 

Lawrencs W. Ga?r, A®?-' » "•2i.-'‘6*^'S’s ^-sbof.-so;*/ 

AVHAOGOM Bessarch sf'd Technoiogy Isbofstorjss 
Anss? Rsi'sarch CsJtter.. iVfoffsU Field, C«2-i>^<?s'fiia 




fe ^y‘ t. / A'-? -'^iv/ 'JyA 


AV^ 


j: 

Nationai Aeronauticr. 

Space A’jministrr.t'ori 

Ames Rasssreti Crfi^ 

Moffett Feld. Caiifc, riis 94n.A5 


j--‘»5^^tff!bH s TAii^gSat A ~~t 

* |K3ts4?s'-^"Ofl UnliB^w 


United Stales Anr,; f' " 
Avicition Rp&ssffCh ^ry f ;0 yj '? 
Devcfoixnant Car'.ta.^o . /f , J 

St. LcAiii,. Missouri ^,31Gu ' 


.1 







ii»my imm Mitif r&M, 6sw:im 

; '>'r ■’’; ^e iimiffri^^ 

fel' ' .V' ■ . 

Taifcar Csl>ge'l'® ind L^jf-siict M. C§n* ^' ■■ 


suifmy 

Rsesrstly a nm procedwr®* wlildi solves the gm^eralng be«ndsry-layer 
equations v/lth Keller's box ngthod, has been developed for* calculating 
"unsteady laminar flows with flow reversal [1]. In regions where the stream- 
wise velocity contains flow reversal^ the solution scheni® was modified by a 
procedure which accounted for the downstream Influence. With this modifica- 
tion, the unsteady flow over a circular cylinder started Impulsively from 
rest was successfully calculated to values of time and space greater than 
In any previous solutions. An examination of unsteady separation for laminar 
fiwi was made and revealed that the unsteady boundary layer for that flow, 
even at large timss, was free of singularities. 

In this report vjs extend the method of ref. [Ij to turbulent boundary 
layers- with f1®-4 reversal. Using the sigebralc eddy riscssity formulation of 
Cefceci and Smith |2]» we consider several test cases to Invgstigate tha 
proposition that unsteady turbulent boundaif^y layers also remain fre® of 
singularities. 


Since tha solution ©f turbulgnt toundsry layers requires a closure 
assurptlon for th® Reynolds shear-stress tsra and the aceursey ©f the 
depend on this ess'OEptlon, we also psrf©« turbulent flw cilcalstlons by 
usiRf th® turbulsnc© rfSdel of irgdsh«« Ferriss and Atell |3j”, w@ solve the 
govtming equations for both madsls by using th@ s©?s® nurerlcal sches® snd 
cospar® the predictions with each other, restricting the cespaHsons to cases 
1n which v?s11 shear is pesitiv 
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Is ®f 1»portiRSt 111 & fiiiMitr of la 4sm^€ still a 

biiff®t1sif iri^l pist stydiss. W^^r* tern of tfes w# f^p^lor tortalofis® 
fiocfels 'ffsplldtl^ ssatsa thst wall shear Is pesfrlvn m4 their ■ 

to tiRstgad^ f1o»$ with flow re^wsal is not aasy. It r^g^ms ^dlftestltss 
to the fuactloaal fonn of the law of the «a11 and to ths in i#t1ch t!i§ ksII 

shear Is detennined. Tvifo n®sr-wal1 assumptions are corssidgred here. In the 
first* the near wall grid point Is located In the logarithmic region and the 
law of the v«ll Is used to link the f1c»/ properties at this grid point to tlie 
walL In the second, a Van Driest formulation due to Cebsci and Smith [21 is 
usecij this Implies that the grid point closest to the wall will occur In the vis- 
cous subl^-'er. 


A further aspect of these flows of current interest is the possibility of 
a singularity occurring In the reversed-flow region. Examples of this phen» 
®venon have also been reported In laminar flows but* In earlier studies febeci 
fl*4j and Bradshaw [Sj have shovi^ that tJj© occurr-ence is not <2 fmtum of the 
governing acfuatior=s but is dus to tte 'liaitatlofiS of tha nissfrlcis,! pi'-ijcedam 
used, k-0 shall demonstrats that* for t?i« osssiples m stu?§y, t^-ssr© ir* no Indica- 
tion cf siidi g singularity In tuj’tulent flow either but there Is & clear Indi- 
cstloFi of rapid thiskenlng of tte bosjndary layer. 

addition t© th® «sf v?s11 functi&-?iSs, t?® to® also soasldtr^d 

two turt'uignc® ssdals for wnstea# fiffitfs without floi^ s'avirsal, TIs® eipbrsfc 
sddy-‘ viscosity forKtulatloi of Csbgsl end Sa^ith (CS| is esf^ared with tli® trans- 
port ffisdsl ©f 8r®dsli»#3 Ferrlss md Atw©11 [3j (BF). Calculations mm p^rforead 
to dstemirss whether ths r®pr®s@ntat1ra ©f unsteady f1©i« with strong prefiSiira 
gradients s-'iqulres that acceunt bs taken ©f transport, of turbyleaci ejusstltfes. 

As will be shown® the pradictic-ns with both fj^dals ar@ stsarly idg-rstlscl for 
both ittady and unsteady flows with and witteut strong prsssur® §ridie?it. 

Ths report has been pr-apsred with six sain sectios'js descrfblna* 
respectively, tha foverning equetlensj the nu!;?er1ca1 procedur®, tha results, 
concluding rer«ar'ks» reftrenoss arsd tha cosjputsr pregrat whleh yse-s enly the 
CS 
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and m r@c8l1 that u' and v' dsnote fluctuations gboyt, the enssK&ls- 
sversgs valocityj u' and v' am zero In unsteady laminar flo ^4 and 
V au/ay Is nsgllgible outside tSia viscous sublayer in a turbulent flew 
These equations are subject to the usual boundsry eor»d1t1on$c ^?1{ich In 
case of boundary layers are 

y 0, u ® V •* 0; y 6 y 


f tha Htynclds strast term, -uV^ introduces m tdditi ?Kia1 

t© ths system §1v@n by Eq$. (2) to {4). In this fn$ present 

calculations using t»?o dlffers-nt tuiisiulsned i©d@1s. fee Is an slgabraic 
eddy- viscosity femii'istlon 4sv@1§psd (for steady floors) fey Cabsd snd Ssith 
e«d tte otter Is & trsrssporl»®£|sist1en ilswlcped by g»/’efi':4i£v«s Ferrlss 
m4 In til® €S 8»ds'l, fe-@ wita Eq, C3) as 
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with two ssparst® fomilss for s^. In ge»ca11ed im&r fmjim% 
bsundify lai^'er C^gj)! dtf!r@d by tte following f©ri?j1s; 
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,li ttli 0 fJlSf f^fteB ‘Is flsffe «4 tilt f# 1 lc«'lsi| firsijis 


s; ‘1 ® I (y •«■■ islsttf 


® b®y?!dsry tetests ths ' <»w?* if»® outer U ©stfMIsSitd fey tfse eestlii” 

ty of the 6 ds|y-*vlsc@s 1 ty forwMlss. 


In the BF R'iodal, wfild-i Is used «i 1 y outside the vIscchis st&leyer, ws 

essufiii T « -sFT^ and v^Hts s sl^gl© first-order partlal-dlffersntlal ©qua- 
tion for it; th® equation was originally developed frosi the tu?i>ul®nt energy 
equation but can be equally wall regarded as an ®p 1 r 1 ca 1 closure of the ®Kact 
shear-strass transport equation. Hil® reads 

3/2 

^ B H- f « U 4- V CtV^) ~ Z&^ (g) 

Her® is a dljiiarsalcnlass auantftye V fs a velocity and L is the dis- 
$1 pet Ion length parisssten specrf 1 «d el^efe«i 1 ca 11 y by l/s » f(n) ^*ith 
n * y /6 end f(n) glv^n frsci an anslytic fit to an sspirfcal cyrvo by 
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n « 0 .. 1 S 

I® < S| < 1.1 


a raore advanced vsf'Slms sf this turhulma aa^l |©| L itself Is «fet®r- 
i®d fm» a transport eq^netioa. 

The turhulent trsisporc velocity ¥^* n^nally (pW ♦ 

prepsrtienil t© a v®1»c1ty seal® sf the large eddies ssiid Is choseti t© b® 
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In am SF a§«li1 «?giiatf<m8, 1fa»r b^<fe<^--cs?M«tims for (1), |S) m4 Cf| 

sr® spplit'd &«fesl4i tfei vfsesi's 's«Ji1«>W9 at y, « S§^/si . In fete 

$tsi#**fl©w sty# repsrtsd Isi |?|s tft®s@ fesujsdsr^f emidft’lsns &mi 
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Hsr® Is svaluated frm the continuity equation (1), end a® Is evaluated 
fro5!t {1} »d (2) on th® assttoptlon that the velocity « Is given by 



(IS) 


f©r 0 < y < y.jj Eq, (13) Is* of e©«ts®» a special 'cass of CIS).-? The ©va1t«a“ 
tl® of Is d1sc!jss@d In Ref. C?|j the lest tern In (IS) css'* be ss largs 
8S hslf tis© sC'Csnd Cpressure^gredftnt) tsvis. In wnst&s# f'iw without fl® 
irevsrsils, ys® the sass liiRer **bsynd®ry“ ccndftiens ®t y^ ® 50 v/u^b but 
tncfmfM of til© pr@s«c® ©f the tlss^fepsndeat tern Ir (2)s a* fet'cesas ®§r@ 
sespllcstsii, If m egs'in Simms that (IS) lislds - reipesb-ar that turhulenee 
stryctur® cf the Infpsr Isy©? Is isnllksly to be sffected unless the extsrnal- 
streffl fRsquancy Is very high - then (1) aid (2) give 
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Kfe «« Keller's tw0“p3f?!!t fin1te-^ffff#»feee fssthod (called tfe Bex 
lotted! ta sol¥s th« syst^a ?*f ®su8t1sr5s ffescrlfe-fd 1« tte pmvims zectim, 
'Ht ®pplle®t1c») of this «®thod te aestesj^ fleas ulth m nm reversal «sir<§ 
the-CS 9®«fe1 has hes^ ctescHted in Mf, |3|, its applfcctiwi to steady twa- 
dlfasmsIciMl fle^s as1fi§ the 8? is dsscHfesd In fef. [?|, Hsm m 
pr®8@nt a Ascription of Urn of ths £S m4&1 to wsteacly two- 

dimensional turtulent flov/s yfl^ t flew reversal as well as a Ascription of 
the ®Ktsns1on of the BF model to unsteady' turbulent flc^i'S with no fltjw 
reversal . 


3.1 CS f4ethOG with and without Flow Reversal 
As In previous studies (see* for exafflple C8J), we transfora the equations 


X = x/Li 


t “ tu„/L, 

0 


n » (u-/vx)’V 


i) 


a dlnsanslcnless stra«» function fCxstiet^t where 


ISli) 


Her© is a refemice ve'Issitys L a r^-ferenc© lenfthj, arisi $ Is t!is 
usual definition ©f the straars fursctlen corr&s ponding to tlj@ ccntirsultj equation 
(1|. With the ralatlons dafinad by (19) and -with the Aflrdticsn of 
viscosity, equations (1) to (3) m4 the boMRdsry conditlens can bs wHtm 


n ® 0, f = f’ ^ 

Prtes denote differentiation with respect to n and 
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(20) 

a*x s¥ 
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f • v'“o 


(21) 




b - H- 


% *“c 

^ \ ® St 


aS w-w 

R! \i 


e 


m 


For slfifsliclty, we shall new drop the bars on x and t, 
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Ms ys® tm «®f>arat 0 selytlsn prefi@dums t® solva tfes $'j%tm glwn by 
&|s* C2S| esd (t1), llfeess tlmm Is m flow Wfvot^al across tft® Isytr* w® ' . 

urn ths staosterd Box. 6a ^ ©timr baad* Mten ther® Is fl^ m^ml# tters '• 
w« MS® thB so-called sif-sfif Bsk t» dsscrffeid 

T© E^s, ftG| and C2‘i) by th® st&BdaH Bsn nsfkmii first »it® 

E^. (20) If# tents ©f thm® first-srder e«jai§t1s»s by intrsKlaclRg fses? dsp®iid«- 
®Rt vgflebtes «Cx 5 »ij,t)s vCxs,nst)s that ISf 

f‘ ** u (23a) 

u ' « V i ( 231 )) 





(23c) 


\4e next consider the net cubs shewn In Fig. 1 and write difference approxi- 

listlons to Eqs. (23). Equations (23a .b) are approximated using centered d1f~ 

fe!*ence quotients and averaged about the midpoint (x^j t„s i.). The 

• n 

difference quotients which are to approximate (23c) ar® ^f^itten about the 
midpoint whose mesh widths as^e 

find hj. This procedure ylslsls the follcxdRg eq«at1cfi:«; 


f. 


- !>« 


.1 






^Ts 


cbv ) j :” ^ 
A«L ^ X 


■j 


X 




^ ^■ f„i * 


' 1 -Js 


« " t , 


n ' n-T 
a 


" b“]'S * ® 

(24a) 

- bit ” " 

(24b) 
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Fig. 1 ftet eybe for standard Box !?ethod. 


ttere hy 


Kan V 'f ths mn of the Vsisluss 

«7 


of V- thrae of 'the fi^ur co-mtrs of the ©f th^ km. Also 

i3 
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1 ... 


*1 " f 




J,n ^ 1 .5. 


The resulting algefcralc s^stes glvsn by Eq*. (24) togetter ^sith the 
boundary conditions » ¥ihich n« becoss 

r “= u« 0, li, “ u,, (25) 

00 a e 

are nonllneijr, W® use totori's tsethod to llr.sartze ths system and solve tbs 
linear system by the block elimination R«thod discussed In rsf, (9J. 
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|£;'-i.sb| fSslI still €®?eW It C?1|S> 

far tis® ioi ©gsitsmd ©t s? 

ppi¥l§as1y mipytsd of I 

Bsk isfthod: If ^rlt© (23e) fsr th® E©k e©»t>sre^ 

Csea Fig. 2) ysing qusntltlss cesttred et P* Q» mi Rj, wfs@r® 




P = Via’ ^ = ^^1“%’ '^n* 


R = Vr 


C£S) 


Ecsuatloo (23s) sen thgn ba wHttsn as 
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y)' (P) •!■ j (fv)(P) ® x( 
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,, .,) * ♦«(!!) II (IS) - ev(Q) II (Q) 


C^) l"^ (P)l 


(27) 
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Fig. 2 Finite difference iriolecule for tbs Zig-Zsg differencing. 





Tf® resMltlRf fs stairs mmUm&r Bt\4 fts Is 

fey astiig pfiCfAi'fs Is tiJ® staBdirtl %m. 


3»2 SF Ifethed wltli m Fits® ilwersal 


Ths scluticn of the gevemlng squitlons for ynstsady f1m% with the BF 
B«5de1, even with no fim reversal across the boundary 'layer. 1s.r»ych «ore dif- 
ficult than with the CS mocteL This 1s because of the hyperbclie nature of 
the govereing equations s together with th® fssniirssar boundary conditions » which 
play an Important role In the solution procedure. As Is cesKKsn in most (If not 
all) mattiods that uss boundary conditions away from the “walls*' the wall shear 
stress Is also an unknown pararster; It can be treated as m gipnvalu® or ss'a 
rachul as described in Ref. [7]. The latter procedure 1s raich P»re efficient 
than the foKisr proeedur® and is used hsf©. 



To selv© the EF iwdsl equatlsns* w® first- Introduce the 


strcKii function 


^(x^y) as In Ref. |7J In order t® sstisfy the sont'inwity ’«1th 

s/t“ 5 w treated as raschult tba resyltlng system can te written as a systea of 
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Again th® system fivert by Eqs. (3a) Is noallnsar and is IlnearlEed by 
ysing toton's siM^thod. This procsdyre gives rise t© the fol Imping fmn 

(2 < j < a) 

® (fg)^ . C2 

”> - / Cfiuj + {2 
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Her© for cmvsnlesse i#s fi8¥© drspptfl tte $ap©r§crtpts Isfl ftav® fefitssd 

CSj,)j (k - 1 . 2 g «.*8 6 ) 


1 ®n 4 




1 - 1 , a. 


Sg)j ® - ®^«j - s^/hj ii-y^ - ^j,4' ^ 







Th© t®rf.@ dsnotad b’r (r^)j (I: ® 1» 2, 3e 4} srs d®fl»isd by: 
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h"h ■’■ 4“% 


For j ® 1» ^ usa tlie 
(15) snd first write thesi as; 


Eqs. (13), (M) 
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u-| ® W|(2.5 1 r 1" 5.2) 
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Aftftr Ks 
eKpresslOTSs wffi get 
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1 .. _ . 123 ^, . ,t ^].1234 


2? “ Cl + j ^1®1^5 '^ "" 2 ^1®n 

y,«-i 

(r^)^ - w^[2.S In -^4. 5.2] -y^ 

,1234. , . » J234.„ c \ 

c.^) "" W 1^1 ““ E^i ) 

u, 1234 


, =* y-j* 


(rj), . («?)”’ - - 4 (? 3 )?;^i ‘ E, - iy,.,EP 
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- Ct^ ~w^ “f ^'s “■■■'? 

F©p J ® J» m MS« the sistjsl bosindary condition » 


In Its 


y, » u 

y e 

iEjy fsm Is 

- Ci^)^ * 0 

The equations (31 ) fs-r 2 < J -s J and 
(32) snd (33) f«i«i a linear system wtilch 1s solved 
@llBlsiatl©n SfSthod dlsit^sssil In Rgf. (§|. 
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y|34) .(i«j 


where <s Is g* pasitiv® csr^stet, ^sleall^ 0lstHfe?st1®a'&?$si?®mt1f’^{'f tV.' 

used by Cefaeci [11 for laminar order to 'study the Cffisputkicn of ■ ■■• <i' 

unsteady laminar flows with flow rwers@t"us1n0 the solution proceslura descHbesI ■■"‘l -c 
In the previous section and to sea vnhsthsr there is a singularity associated^ ■ ‘ -' 
with such flows. 


In performing calculstlons- for this case and for the ©tters considered. -■ 
here* cars must be taken In gangrating the Initial condltlws 1n tte (t^y) '■ 

and Cx^y) planes at so«a distance*- ssy % « For a laalnar flow If 

® Os tha initial velocity profile -for the velocity distrlbvttsfs given fey 
£q. (34) can.be taken as Blaslys 3*sd there' 1$ -im? difficulty ateri ca*p«jt-if5§ " 
the sslutfon 'In js > 0 since th® Initial bowdafy layer Is ©f zero th1d®ess. 
If n^fQt m can take 

« 1 “> £s{Kp “ Kg)(t^ t^) ^ ■ Q < K < 

but tten ®t K “ K^- thfre Is & d1sc®f.t1nu1ty In ttes prassur® ff-MIsnt, Since 

it acts on alrsady-estebli^td- bo^cfev^ ‘ Isyer* Initial ysspsas® is 

inviscid leading for^lly to ® v©1©s1ty slip mi a sttocun^i^ l^yev* at 
the wall. 'The' treatment of tlit bouadlary layer is tf?§n rather ss&tle fsee -Ref.- 
IIOJ) but if we are not too eofmtm4 with the details of' the'SolBtieB m&r 

jc ® Kq 9 which is th® cesa hsr®* a cs^^^anlirst pmc©dist^i.fe®i1d b© to writ® 

Eq. (34) as ''; ' 

L-v . ■■ 

» 1 -<• aF[(.x - KQ)/i.Hs - s?Ht^ - t^5 (35) 

whtre F Is a sr:;0Oth fsmctlcft «h1di'vs,s1shas If js < snd Is «?s1ty if 
K — K > a. For exfiiTipl®, ws can taks f(s) ® slnCws/S) ©■'< s < I* sad 

0 ‘ f 

a » 0.05 with ten stations t^t^s«a -x'-* ». and x » x,"+'a. A sisiilar dif- 

fleylty would ©cc^r at t ® 0 If. ‘in Eq. (34) v;sr@ replaced fey t 'sine© 
th® boundary layer Js well ©stabllshed it t » 0. 
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As can b® sgeos beth ®Kpr®ss1e??s give nearly tha ses® results. 



Th® results In Fig. 4, as In laminar flo-rfs, gKhiblt no signs of singular- 
ity for all calculated values of t. This is In contrast to thB findings of 
Pats! and Hash [11]. Again, as in laminar flov/s [3], we see the familiar 
rapid thickening of the boundary layer In the reversed flo*#? rsglon. If it 
had not been for this, the calculations would have been ecssputed for greater 
values of t than those considered here. 


Tha two other test cases eonsidared here correspord to C??s«s 4 and 5, 
as reported by Carr [12]. Case 4 Is fcr unsteady Howsrth flew, It starts 
from a wall “established staiidy flat-piste flc^* cn whic« s 
tion of Ug is Imposed at t ® 0. The external velocity dlstHfcistlftti is 
given by 

Ug “ 1 ~7(m ~ 1,24)t 1.24 £ St < 4,m (37) 

' ».i ,.1 

wters. e Is a constant equal to 2 A/3,45 sec m Th® fisii. was assus'sid tc 
b@ steady up to k ® 1.24r.p th® velocity dlstrlbytien Eq. (37) «as tf?gn 
Imposed as ® furiction of s adid t. This test case d1ff@i:s fre-a th® p’mvloivs 
on® In that, once ths flow soparstss, it dees mt reattsdi. For this reasosi, 
the csicylations' can only be continued as far as th® ststitn whsro the flew 
rovers®! firet occurs. Tte initial velocity profil® at s *® 1.24- and for 
all tisiis correspond to a f1at"plat® profile with a wmentuis ihir.krsess Reynolds 
n5S!b@f (Rg) of 4860, and local skIn-fHction ccefflcient C|, of 2.8 x 10“"\ 

As In ti» .previous t®st casst, we Introduce a functira F, so that at 

5 

X ® 1.24, dUg/dx ® 0. Since w@ a1§o want the solutiens t « 3 to cor- 
respond to ste?4y-stat@ solutions, m iRtroduc® ®sther fy^K-tlors Fg ordsr 
to Sit 8«g/st ® 0. mth ttes® ftsietie-?js, £c|. (37) tfesn bec^sss 
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^ ,.«r I 9 -a, 

S & %S>! f 


® slf3 


Pipm 5, $ md 7 nk¥j tit 6f,1e»1tt®«l 1®et1 sMn-fHctlcs^ emfflcimt 
C|s tte tii&pf f2«6«r 'N te# tl?© psfitstssi llikisifsg P^yitslis mMwf' ■ 7^" fsr 
iHH ttst ces®. %f es1ol6tie«s ^sni d&is by niliif tefcfi CS m4 SF B©d@1si 


ths msislts sho3vis by solid '\im% refer ft© tte predictions ef tte CS s©ds1 and 
t«»0S8 shWR by circles refer te the predict'iefjs of tli© SF w4s1. 


As seen from thasa threa flgijreSs there Is essentially no difference 
bafe^ssn the predictions of both snodsls. Altiwush them Is soss discrepancy 
In the shgpa factor pmd^ctiwis, this dots mt s®©i to b® tea significant. 






Figure 6. Ceapyted sh»pa factor distHbutlosj 


for test CuS '2 



ne«r« ?» CffiRputsd iBKi«ntss»tl?felxess BsfBsIds msMr i%t t,?st s;is® 4.. 






" AeeoNIng tko pmdietims ef tl'« €S ®1sa itm tte cKpafelV 

itf ef pmlfctiiif bmnd&ry ley&rs wi^ flm rwml* fefe® 

¥»lstes first srs®4 t O.ES, x » 4.®9* Si&® tfe§ §f M-adarf 

IC^trg f@r %*a1«#s sf k la tli© nesgs 1.24 ^ k < i,®§ f©? t > C.fl tosstls 
Wss sptelflestlea ®f s wricsfty preffle at a » 4.SS, f^®nta smh & 
preflle Isy assslrtg It Is giwn isy tfsa estra^1at1@-?8 ©f ta® t^l«sclty prefll^s 
for x < 4,S3c Thls j*n?et®ym ie t^idi tte gstrgps1ate<^ statlssi 
M a d«nstr®aa boanctery c®!6lt.fs»is. s1l«$ ths ealculitloE^s t® te cestljss^sd l?i 
tte negative mil shear region as slm'OT In Fig. 8. 


The third case consit^red In our study corrasportds to Case 5 irt ref. [12], 
which If) a way resembles the estsrnal velocity distribution In Eq. (34). It 
Is given by 

Ug * 1 + {Btr [C-Cgft ■»• (39) 

-'1 

where A » 0.05, B = 3.4 sec , 5 » (x — 1.24)/3.45 md the range of x 
values are limited to 1.24 £ x < 4.69. As bsfos'^, the Initial velocity pro- 
files at X = 1.24 for ^11 t correspond to a steady flat-plate flow with 
Rg « 4860, ® 2.89 x 10”"".. W« again rsodify £q, (33) to avoid the dis- 

c^itlfiulty in tfte pressure gredlest. 11)1$ tifss we ?jiult1p'?y tte right-hand 
side of Eq. (39) by F-j used I 0 £q, (38), 



Figure 8, ^tristim of wall-shsar pamsster f” with distance as a function 
of tte for test ease 4. 
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fim^m f 10 ilrai Urn «s1e;JakisJ 

mi ^ f®r tills tet ««;■?'. 

m ^m&mt tte pjf4l€tl@is 4vf fe^ti mists.- f’lpr«.ll fev-xis 

vsl@sltf. preflles- fer s iiii ' E»sttSf®iif* ar.a 

fr» teess fl®«p@s, tM'pri&4letlw &f astfe ttirlMitei« i« fis? s$m 

far 


ill® wHatlffii ©f «sa11 sli€®r §®rwst«'r f** ®g s 
ef « ®?id ts ifid Flpre 13 sh»s tte «1§sfty prefllis® fee1i,i^» 

Ing tha reglats In vi’.ldi tt.2P8 Is flea pav«rsa1 moss tendary laf&r. 
coTifJutetlcns which ar© dons by ysing th@ CS srsdsl prml^ cmfirmtim 


of the ganeral trends In test esse 4* nm^ly that as in laalnar flsa-Ss 
the mstB&d^ turbulent boundary layers ttilcten rapidly isitf! lRCt‘'@as 1 n 9 flw 
reversal. A nw feature h«@ver Is the dip in tte grgphs ef f® reiar 
s » 2.§ which develops as t Increases towards 0.40. It Is pessibis tint a 
singularity occurs In the solution at a later tfcs ®s fiany Mihers liev@ su§™ 
ftsted Is the case for laminar beundary layers. Tit® mmt eoptit arg!.s«tit In 
favor of this ^jenaisnosTi has bam advanced by Shen flSf but m nota tbst 


tm Rjost daiinit© sign sf its ©csirrK/ias eppsared In his grapfe ®f dlip1nc4®sfit 

thickness wbldi showed spikey characteristics.. Hero tte dispisossant thickness 

s®€8is to be fairly saicoth but the skin fristica bea»s splfcsy. 
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Figara 9. Ccs?i> 5 ited loal sJ;i^-fr'1ct1c«n diitv^ibutian for tsst cass S. 
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Figure 10. Craputed i?»;»:fitffi!“thickn®ss Rayriolds nuinaer for tsst case 5. 
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Flpr^ 11. Crasparlson of calculated velocity profiles for test case 5 
with rsa f'lw reversal. 
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FI gyre 13. 


Calcislwtsd velcslty profiles Includlss fl®? reversal by ths 
CS Bwiels for test case 5. 
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f, ctiaiMM mmm . ' ' . 

Based on tha stadlts c^s^selsd 1» tl?l§ rsp§i*ts i*®'&teswt that: 

1, Tift mifirlsal sal^tlss ©f «mst©My iMfw ®M t«flja1f?5'it fe«jr^«ry 
layers Including tli® f1®» revorstl ssrsss tM laytf s« fes slstalr®« 
quits MtHfectcHly for s ff¥t^ pi^ssum A' urina- 

tion cf &stli rT'ful^^p 8Ri zlg» 2 eg tox sclss^s sra sl>®w?i te yield 
accurate results for urssteg# fess^ndary layers. 


I? 


Z. Mhethsr the unsteady bouiidary- layer ©qyatlons fsr laminar and turbu- 
Isnt flc« are singular for a givsn pressure dlstributlor: still 
remains to b® invgstipted. TIss results for test cas® S incllcat® 
that at large tisss ther® is a puKling "Id^k'’ 1n the wall stesr 
parareter> this i»y be du® to a sis^gulartty or it sisy b« due 
to a ?si«rlcs1 problesa. studies conducted by CctecI |14| and 

van OsKsslen and Sh« f15] for & eircwigr cyllfidtr started 
iMpulsIvely fross rest Indicate that ®t lar^© tlsfisss t * 1.25 or 
Mr@® tlisre appears to be a sifigylsHty In ar®«d ^ ®* 120®. 
M©^@¥@f 9 th®s@ ctlai1atio:ss d© not iRdlcsts sij’sy pux'jllfig b&hmiw 
in thttal'i $hmr psrsiMter near “slafylarltyj” tht f“-vatues 
srostli md ^11 bthavsd fsr tt©ss and liirftr t lists, fts tbs 
(ftWiv hfip.ds «m 1 f 5 lrrg tli® s'«'»r®s«1ts fer t@st eis© S» t« find tl'iat 
iilill® thm® Is &n a&r?§rsa1 teharior In it lar^s tteSs tli@ 
cofpsspsadi ftf S^“%*a1uss ®re ^©©th sad fe^sll behi¥@d» ® treiid Khldi 
is QBs^site t@ that far i el real sr cylls^dar. 


3, A cffispsHsen of tte prsdietlens of Vm tuffeu1e?sca ®d®1s» ftarsly* 
CS md SF Indicate that for attacliid f1©iSs, both mdsH 

yl@1d alifsst Islsritlcsl results. This Is also tr»j® fer fle^s 
»h1ch are s«ff1deat1y stnsrjg In pressure §rad1ent t© catis© flow 
rsvtfsal across tte layer. 
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®?#mt!sTif Ifea tepi^t- 1@ ttes c^^spiter pf^ffi^i ©assists sf f»Mf tFPss »f 

.csF^, Cari 1 €®Btsls‘$§ tlss felt!® ©f tfea f 1 e» psft^ls® s®dsf ©mst^i 
ten! t fm0ms l&s faTififtssf Irafasmti® 
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1 itystmr @f t“Stat 1 ©ss ts bt ea 1 cy 1 ®tgd 
©f K“Stst 1 em te ba calealatsd 


fITR ^“Ststlan where transition begins. If ths Inltlel velocity 

profile is for turijulent fisiafs* than HTR« 1 . If fim is all 
1 m 1 nar> set EfTMZl. 

IBD¥ Sp-scifiss v^hatter the flo^ at k “0 starts as a flat-plate 

flK'^ or as a sta§nat 1 eaj-po 1 nt flo^, 

®1 f 1 at“p'(ate flew 
°2 stagnation-point flow 
RL Frse-strsafii ^ynslds ni^er® u^L/v« 

IPfjff Controls the print iratpyt 

«1 prints ©ut only ths feow» 4 ary- 1 sy£r paresst'srs 3 *>o, c~« 

v 

Rgj H m4 external valosfly dlstrlfeutlo-n. 

» 2 ' prints cut profiles as «»11 ts the botjndsfy-leytr 


DSTACl) &M ¥SP are the ncr^silfom §Hd 
teress t^is layer, Ift® grid used In this report 

t the rati® of 1 ®g^s ©f any 


censtanti that ISg &nj ® distance 


rs tSist control 
® mmMric pro§r@ 
adjacent fnt®' 
the j»th line Is ^ 



the foil owing 




There are tw® para®ters In this t?^u£t 1 en: » ths length of tte first steps 

Md Ks tis® ratio ©f twQ successive steps. The tetsl nu^Jjsr of points J 
eati fee calculated frea the following fomuls: 

Inp 4. CK- 1 )Cy h, )] 
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K art dWiS-gfi witfe tj-ffcal fm' m4§rsm ®f ail ®«i 

1.3sf Is fiSitrals ippr'snig^tslf f§ .|pld tf-a fe-.sdar-j' 

l®l?«r tr®. siiffl€fs«t t§ %mAm? 

SiSMii mlm& @f li| mi I fensst fe§ s?#jli Urn 
m meimr e§ rli ssSst msm4im ts « fltiti §? g-si^iiitc- 3 7?^. '£|„- 1 | 

Mt §11ea J t© exsisd 1Q1. Fifyrs 14 Is prmi6&4^ tlisrf f@rs * t© pwlds 
gaftec® In thg sslectl&a ©f 4. 



Flpre 1#. ¥srfit1c?! of K t^ltfi for n -faltiis. 
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CF and ilTK &m the local wd values %McJj am t» start the 
turbulent f1®.« e§1cy1at1©ss» Ths liiltliil velocity sjmffi® Is gtsisrated by 
usfef the foifwjigs proposed fay Sresivllle fsis ref, S) 

. j fl it’^j 

vr* SB ^ |in y c nC^ cssmi (n'^ “** ??“ll (4S) 
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frm Eq. (40| and fr» tte daflnltlcns of 5* and §» It «:&« be sh«p thst 
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K and then substitute then Into Eqs, (40) and (41 )s thus 

«“Profi1s. 


Cards 3 and 4 r&acJ m the t end k stoticr.s, mspscttvsiy* TT» 

pressiit cKsputer pnj’ffs^' s^icifies thi sMtarnal ViSlccIty dfstribiitioa &y s 
f®mi1ft and coaputss the .^in^ensi^siless pressure grad1e»it parasasters s?m1ytlcffi11y 
as is sSm® in ths listing which f@11®s this ssctleh. Tte test case in the 
coTfuter listing Is fi^r c^^s® 4 §f ref, |12|, 


Ojitpiit 

?k'ps?5dln§ Qff tte IPfJ'IT, tte esspytar pr§fria prints 9«t ths pmffles 

f. f‘# f* fflid b as a fyfiCtI® of the sImfleHty vsHafele n and grid 
parssster j to^thr-f* with a psn»t®r KM.C(«I). Iter®- lfMJC0) « 0 Khan 

us use tte stgndani and ® 1 vii’sen v?® sass tte slg-zaf bsss. 
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